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The low-temperature coarsening dynamics of a one-dimensional Ising model, with conserved 
magnetisation and subject to a small external driving force, is studied analytically in the limit 
where the volume fraction /j, of the minority phase is small, and numerically for general fi. The 
mean domain size L(t) grows as t 1 ^ 2 in all cases, and the domain-size distribution for domains of 
one sign is very well described by the form Pi(l) oc (l/L 3 ) exp[— A(/i)(7 2 /L 2 )], which is exact for 
small /i (and possibly for all fi). The persistence exponent for the minority phase has the value 3/2 
for fi — > 0. 
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I. INTRODUCTION 

The field of phase-ordering dynamics is by now quite 
well developed [fij. It deals with the approach to equi- 
librium of a system quenched from a homogeneous high- 
temperature phase into a two-phase region. Familiar ex- 
amples are binary alloys and binary liquids, which are de- 
scribed by a scalar order parameter. Recent work has ad- 
dressed also cases where the order parameter symmetry 
is continuous rather than discrete Q . Especially interest- 
ing is the scaling regime that emerges in the late stages 
of growth. For a scalar order parameter, for example, 
the domain morphology is apparently time-independent 
if lengths are scaled to a single time-dependent length 
scale L(t), which represents the typical 'domain size'. 
This implies that two-point correlation functions depend 
on the spatial separation r of the points only through the 
ratio \r\/L(t). 

By contrast, the coarsening dynamics of driven systems 
has been much less studied. A physically relevant exam- 
ple, indeed the motivation for the present work, is the 
phase separation of a binary liquid under gravity. Nu- 
merical simulations of a (physically less realistic) alloy 
model with gravity suggest the existence of two growing 
length scales, parallel and perpendicular to the field g] 
- pi], although it has proved difficult to unambiguously 
extract the time-dependence of these length scales. 

An independent field of study concerns the stationary 
properties of these 'driven diffusive systems' ||. Here we 
focus on the nonstationary, coarsening dynamics which, 
as we have noted, has attracted relatively little attention 
thus far. Ultimately we would like to understand the 
coarsening of binary liquids under weak gravity, includ- 
ing hydrodynamic effects |(| , but as a first step we settle 
here for a less ambitious goal. Specifically we study a one- 
dimensional Ising model with conserved dynamics and a 
driving field E which favours transport of 'up' spins to 
the right (and 'down' spins to the left). We work in the 
regime T <C E <C J where T and J are the temperature 
and exchange coupling respectively. We derive exact re- 



sults in the limit where one phase occupies a small volume 
fraction, and numerical results for general volume frac- 
tions. The main results are a yt dependence for the mean 
domain size, and a domain-size distribution for domains 
of one sign of the form Pi (I) oc (l/L 3 ) exp(— \l 2 /L 2 ), 
where L is the mean size of domains of that sign. We 
also show that the recently introduced 'persistence 
exponent' 6, which describes the fraction f(t) ~ t~ e of 
spins of one phase which have not flipped up to time t (in 
a sense to be clarified below) is 9 = 3/2 for the minor- 
ity phase in the limit where that phase has a vanishingly 
small volume fraction. These are the first analytical re- 
sults for the coarsening dynamics of this driven diffusive 
system. 

The paper is organised as follows. In the following 
section we define the model. In section III we discuss do- 
main growth and dynamical scaling in the model, while 
section IV deals with the persistence exponent. Section V 
concludes with a summary and discussion of the results. 



II. THE MODEL 

The microscopic model we consider is a chain of 
Ising spins Si — ±1 with nearest-neighbour coupling 
strength J. The system evolves by nearest-neighbour 
spin-exchange dynamics, with a driving force E that 
favours motion of 'up' spins to the right over motion to 
the left. That is, the microscopic processes are 

++ — ^ +-+- A = 4J-£(i) 

— ++ ^ -+-+ A = 4J + £(ii) 

++-+ ^ +-++ A = —E (hi) 

-+ — r± — +- A = —E (iv), 

where the rate for a process from left to right is propor- 
tional to (l/2)(l-tanh(A/2T)). We distinguish between 
the 'forward' and 'backward' versions of the processes de- 
picted above by using — A and to denote the process 
from left-to-right and right-to-left respectively. 
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We consider the regime T <C E <C J. This is very 
different from other studies, which have concentrated on 
the limit J/T — > ||. The system possesses metastable 
states consisting of long domains of parallel spins, sep- 
arated by domain walls. After a long time of order 
cxp[(4J — E)/T], either one of processes (i— !l ) or (ii— !l ) 
takes place, i.e. a spin splits off from a domain. If process 
(ii — x ) takes place, the system quickly relaxes back to the 
metastable state by the reverse process (u\— ), since all 
the other possible processes are endothermic. If process 
(i— occurs, the system can relax further by exother- 
mic processes of the kind (iv— *■), so that the 'up' spin 
moves to the right, eventually meeting and adhering to 
a domain wall; the system may also relax exothermically 
by processes (iii— *■), so that a 'down' spin moves left, 
eventually meeting another domain wall. The motion of 
these free spins is unidirectional, because the reverse pro- 
cesses (ihV- ) and (iv-*— ) are inhibited by a factor of order 
exp(-£/T). 

The result of the free 'up' spin moving to the right is 
for the 'down' domain through which the spin has trav- 
elled to have moved bodily one step to the left; the result 
of a free 'down' spin moving to the left is for the 'up' 
domain to the left of the wall to take one step to the 
right. We may therefore map the microscopic dynamics 
of the lattice of spins onto one for an array of domains. 
The system then evolves by a domain of 'up' spins mov- 
ing spontaneously to the right, or a domain of 'down' 
spins moving to the left. The rates for such processes are 
independent of the domain size. This mapping is anal- 
ogous to a mapping by Huse and Majumdar M for the 
low-temperature Kawasaki chain, corresponding to this 
model in the opposite limit E = 0. 

When domains are of size two spins or less, they can 
vanish. The microscopic mechanism for this in the Ising 
spin picture is series of events involving two random walk- 
ers that may coalesce, which translates into rather com- 
plex transitions in the domain representation. However, 
when the domains are large, the details of this domain 
annihilation process are not expected to be important, 
so we choose to study a model where the simple domain- 
shifting dynamics applies for domains down to size one, 
and removing a domain if its size reduces to zero. Sim- 
ulations of this simplified system permitted much better 
statistics than would be possible with the true, micro- 
scopic system, whilst still giving indistinguishable scaling 
behaviour. 

The algorithm used for simulation was the following: 
(i) set up an array of alternating 'up' and 'down' do- 
mains; (ii) pick a domain at random; (iii) if the domain 
is 'up', move to the right [i.e., reduce the size of its right 
neighbour by one, and increase its left neighbour by 1], 
otherwise move to the left; (iv) if one of the neighbouring 
domains is of zero size then remove it, merging its neigh- 
bours; (v) update the clock by 1/ (number of domains); 
(vi) repeat steps (ii)-(vi). 



III. DOMAIN GROWTH AND DYNAMICAL 
SCALING 

A. Simulation results 

Since the number of domains in the model can de- 
crease but not increase, the average domain size must 
increase monotonically. In a finite system (with periodic 
boundary conditions) the system will coarsen until there 
is only one 'up' domain and one 'down' domain; this state 
will not be stationary, because the dynamics still permits 
both domain walls to perform correlated random walks. 
There will be a wide regime of time during which the av- 
erage domain size is much smaller than the system size, 
and the system might be expected to display dynamic 
scaling. 

Simulations using the domain model described in the 
previous section were performed, using lattice sizes in the 
region 10 5 -10 6 spins, for times up to 3 x 10 4 and averag- 
ing over several hundred samples. Random initial condi- 
tions, where a spin has a probability /i of being 'up' and 
(1 — fi) of being 'down', were used; similar results were 
found if an ordered initial state was prepared using alter- 
nating single 'up' spins and domains of (1 — 'down' 
spins. The dynamics conserves the magnetisation, so the 
volume fraction /i remains unchanged. Several different 
random number generators were used, and the results 
checked for consistency; the four-register shift generator 
of Ziff jl(| was used for the runs of highest statistics, 
once it had been established that it gave results consis- 
tent with other generators. 
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FIG. 1. The average domain density N(t) for four values 
of the volume fraction fi. The straight lines are plots of the 
asymptotic prediction from Eqn. (lid) 

Figure [j] shows the average domain density plotted, for 
different volume fractions [A, as a function of time t on a 
log-log scale. The straight lines all have gradient —0.50. 
Figure [2] shows a time-dependent effective exponent, de- 
fined as the gradient of a line between successive points 
in Fig. |[ plotted as a function of l/ln(i); the results 
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show that the data appear to approach the value —0.50 
as t — > oo, though the convergence is slower for larger 
values of /i. The characteristic domain size therefore in- 
creases like t 50 . 
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FIG. 2. The effective exponent for the decay of the domain 
wall density. 




FIG. 3. Scaling plot of the domain size distribution for vol- 
ume fraction /i = 0.5. 

Figure || is a scaling plot of the domain size distri- 
bution Pi (I) (defined as the number of domains of size 
I per lattice site), for a simulation with fi — 0.5. The 
average domain size L(t) is defined by L _1 = J^i Pl(J)- 
The data show good collapse to scaling, even for short 
times. Figure || shows the same data plotted in the form 
ln[L 3 Pi(l)/l] versus (l/L) 2 . The linear behaviour evident 
in the plot suggests scaling of the form 



P cx 



exp 



{-\[l/L(t)f} 



(1) 



For the case \i ^ 0.5, the average sizes of 'up' and 'down' 
domains differs by a factor fi/(l — fi), but nevertheless the 
up- and down-domain size distributions P + and P_ are 
both found to satisfy independently scaling of the form 

m. 



10" ; 




t=10 
t=100 
t=1 000 



10 

(l/L) k 



15 



20 



FIG. 4. Fit of the data in Fig. | to the form (Q) 
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FIG. 5. Scaling plot of the equal-time spin-spin correlation 
function for p = 0.5. The prediction assuming uncorrelated 
domains is shown for comparison. 

Figure |^ shows the scaling of the equal-time two-point 
correlation function C(x,t) = (Si(t)Si +x (t)). Although 
the data appear to collapse to a scaling form, the ap- 
proach appears slower than was the case for P;, suggest- 
ing that Pi (I) is a more natural quantity to describe the 
system. The apparent simple form for Pi suggests that 
the scaling state might be very simple, for instance there 
might be no correlations between domains in the scaling 
limit. The structure factor S(q) [the Fourier transform 
of C(x)] for a system consisting of uncorrelated domains 
may be shown to be [O] 



S(q) = 



Lq 2 



(2) 



where Pi(q) is the Fourier transform of Pj. The inverse 
Fourier transform of (||), where P; was calculated by as- 
suming the simple form (|l]), is plotted in Fig. [5]for com- 
parison. The discrepancy with the simulation data shows 
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that strong correlations need to be taken in to account 
in this system, even though simulations measured only 
~ 3-5% correlations between the sizes of neighbouring 
domains. 



B. Solution for domain density assuming scaling 



N(t) 



{Trf Di- 2 + (l-/i)- 2 ]} 



1/2- 



(10) 



The straight lines in Fig. [l] are, in fact, plots of equation 
( |Io| ) for appropriate values of /i. The excellent agreement 
of the data with the prediction confirms both the pre- 
dicted /i-dependence and also the simple form for F(x). 



The numerical data suggest that the scaling function 
for the domain size distribution is the same for minority 
and majority domains, and is independent of fi. Using 
this as an assumption, we can calculate the average do- 
main size. The number of domain walls changes by a 
domain of size one shrinking to nothing, and so we ex- 
pect that P±(l±) ~ l± for l± — * 0. We therefore assume 
a scaling form 



(3) 



where a± is chosen so that F(0) = 1, and F'(0) = 0. 
The rate at which 'up' domains vanish is P+(l), which 
approaches dP+/dl+ in the continuum limit. Since each 
domain vanishing event removes 2 domain walls, the rate 
of decay of the number N (t) of domain walls per site is 



dN 
~dt 



= 2 



( dP+ 




dP- 




\dl+ 


l + =0 


dl- 


i-=o) 



= 2 



L\ + L 3 _ 



(4) 



The domain size distribution is normalized to the den- 
sity of domains, and there are the same number of 'up' 
as 'down' domains, so 



n/2= r P+ {i + )di+=^ 

Jo L + 

roo £ 

N/2 = J P+ (l_)dl_ = -^ 



(5) 
(6) 



where /i = xF(x) dx. The densities of up and down 
spins are 



A 1 = 



I'OC 

/ l+P+(l + )dl + =a + h, 
Jo 

POO 

l-/x= / J_P_(Z_)cM_ = a_/ 2 , 
Jo 



(7) 
(8) 



where fa = J °° x 2 F{x) dx. Substituting for L± and a± 
from (|]-||) into (^), and integrating, we find the following 
asymptotic result as t — > oo: 



N(t) 



/|t[ M - 2 + (l-/i)- 2 ] 



1/2 



(9) 



For the particular case F(x) — exp(— Aa; 2 ), suggested 
by the data, we have /i = 1/(2A), fa = tt 1 /^" 3 / 2 / '4, and 



C. Solution for P± in the limit fj, — > 0. 

The simulations are in excellent agreement with the 
scaling hypothesis, and with the simple form for the scal- 
ing function ([!). We would like to have an ab initio ex- 
planation for these results. Unfortunately, we were only 
able to solve the dynamics in the 'off-critical' limit [i — > 
[or, equivalently, fi — > 1]. 

In the limit (J, — > 0, it is necessary only to consider 
the motion of the minority spins, i.e. the motion of the 
majority domains. This is because the size of each do- 
main performs a random walk until it either dies (shrinks 
to zero size) or coalesces with the nearest domain of the 
same type. To coalesce with another domain, the inter- 
vening domain of the opposite type has to shrink to zero 
size. The time scale for the vanishing of majority do- 
mains will therefore be much longer (by a factor ~ M ) 
than for minority domains. The dynamics will therefore 
progress primarily by minority domains shrinking to zero 
size, and never coalescing, whereas the majority domains' 
sizes only change appreciably due to coalescence. 

Consider a particular minority domain containing n 
'up' spins at time t. In the limit where coalescence is 
forbidden, the domain changes size by: (a) an 'up' spin 
arriving from the next domain to the left; or (b) an 'up' 
spin splitting off the domain, and moving to the next do- 
main on the right. It is also possible for a 'down' spin to 
move from the right to the left of the domain, but this 
does not change the value of n. Notice that the dynam- 
ics is independent of where the neighbouring domains 
arc, whether they are vanishing or coalescing. Once the 
size of the domain reduces to zero, the domain ceases to 
exist. The Master equation for the probability P(n, t) of 
the size being n at time t is 



dP(n, t) 
~dt 



= P(n + l,t) + P(n-l,t)-2P(n,t) (11) 



for n > 1, with P(0,t) = 0. In the continuum limit 
this Master equation approaches the diffusion equation , 
whose solution with P(n, t = 0) = S(n — no) is 



P(n, t) = (Ant) 



-1/2 





(n- n ) 2 


jexp 


At 



exp 



(n + n y 
At 



(irt) 1 sinn exp , 

y ' V 2t J V V At At 



(12) 
(13) 
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In the scaling limit t 
this reduces to 



P(n,t) 




In a random initial state, the number of domains of size 
no per lattice site is (1 — /i) 2 /i"°. The 'up' domain size 
distribution at a time t in the scaling limit is therefore 



p+(i+) =5>-/i)V°P0+,*) 



(15) 
(16) 



which is of the form @, with L = It 1 ! 2 . The total do- 
main density (twice the density of 'up' domains) is 




N(t)= r p + (i + )di + =^= 

JO V7TI 



(17) 



which approaches (|10|) in the limit ll — > 0. 

We may calculate the size distribution of the majority 
domains from the probability that a given region contain 
no domain walls. Consider a region of M lattice sites 
containing n minority spins; these spins need not neces- 
sarily all belong to the same domain. Then the number 
of minority spins changes by spins entering the region at 
the left and coalescing with the leftmost domain wall in 
the region, and spins splitting off the rightmost domain 
wall. Once n reduces to zero, however, any spins enter- 
ing the region from the left will simply pass through the 
region, and it will remain empty. We have again assumed 
the limit fi — > 0, by not allowing for any domains to move 
out of the region or to coalesce with domains outside the 
region. 

Under these conditions, the probability Putin, t) obeys 
a Master equation of the form ([ll]), whose solution with 
n = jiM at t = is 



P M (n,t) = (Ant) 



-1/2 



exp 



exp 



(n - fiMf 
It 

(n + uM) 2 
4i 



(18) 



The probability Pe of the region M being empty is then 



P E (M,t) = 1- p 
Jo 



m (n, t) dn 



(19) 



U 2 \ 



= 1 - (47rf)- x / 2 J ^ exp du (20) 

Consider now a region of m sites, sitting inside a domain 
of size I sites. The number of positions that it can oc- 
cupy within the region is (I — m)Q(l — to), where is the 
Heavyside function. The probability that a randomly- 
chosen interval m lies within a domain of size in the range 



I to l + dl is therefore (l — m)Q(l — m)Pi(l) dl, so the prob- 
ability that a region of size to contains no domain walls 
is 



P E (m) 



(l-m)Pi{l)dl. 



(21) 



Differentiating twice, we have P^(m) — P/(m). 

The result for the distribution of majority domains is 
therefore 



P_(L 



d 2 P E {M) 



dAP 



M=l- 



H 3 L 



: exp 



(22) 
(23) 



which is of the form (|T|) with L_ = L + / ' [i. This shows 
that the calculation is only valid to lowest order in 
since the conservation of magnetisation implies that 
L-/L+ = (1 - ,j)/ M . 



IV. PERSISTENCE EXPONENT 



It was recently found that, in a ID Ising model evolv- 
ing at zero temperature from a disordered state under 
Glauber dynamics, the probability that a given spin has 
never flipped up to time t decays as t~ 6 M. The value 
of the 'persistence exponent' 6 depends upon the mag- 
netisation, and whether the spin is in the majority or 
minority phase. The dynamics for a spin in a phase with 
volume fraction p, is the same as for a q-state Potts model 
with symmetric initial condition, with q = 1/fi p2[ , lead- 
ing to the result that 9 takes values in the range 0-1 as 
ix decreases from 1 to 0, with 6 = 3/8 for ll = 1/2 [ ^3| . 

For a ID driven diffusive system, there are two kinds 
of persistence that may be considered. The first con- 
cerns the probability that a spin in the microscopic Ising 
representation never having changed its value. When- 
ever a domain wall emits a spin, that spin moves rapidly 
through a domain, causing each of the spins in that do- 
main to have flipped twice. Since this spin emission is 
a Markov process, the probability that a given spin has 
never flipped decays exponentially. 

A more interesting kind of persistence to investigate 
is the probability that a given spin has never belonged 
to another stable domain. That is, we discount the 
rapid flipping due to spin motion through a domain, and 
consider only the case where an entire domain has mi- 
grated towards the site in question. This kind of coarse- 
graining in time would also be necessary when studying 
the Glauber-Ising model at low but non-zero tempera- 
tures, where short-lived thermally-activated flipping of 
spins within the interiors of domains occurs, in order to 
recover the true zero-temperature persistence behaviour. 
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A. Analytical results for /i 







Let us consider a test site initially within a minority do- 
main, a distance n% sites from the left domain wall and n 2 
sites from the right wall. We define ri\ and n-i such that 
a spin in a domain of size unity has (ni, n%) — (1, 1). The 
dynamics causes the domain walls to wander stochasti- 
cally, and eventually one of the walls will cross the test 
site, i.e. the test site will have nipped. 

We shall assume that this domain does not coalesce 
with another domain before one of the domain walls 
reaches the test site. This assumption will certainly be 
valid in the limit fi — > 0. Then the three processes that 
cause the domain walls to move are: (i) an 'up' spin 
joins onto the left-hand edge, (ni, 712) — * (tii + 1, 712); (ii) 
an 'up' spin splits off the right-hand edge, (711,712) — > 
(711,712 — 1); (hi) a 'down' spin splits off the right-hand 
edge and moves through the domain to the left-hand 
edge, (711,712) — ► (ii — l,7i2 + !)• The Master equation 
for the joint probability P(n\, n 2 , t) is then 



dP(n 1 ,n 2 ,t) _ 
dt ~ 

which becomes, 

X 2 , 



P(n 1 -l,n 2 ) + P(n 1 + l,n 2 -l) (24) 
+P(n 1 ,n 2 + l)-3P(n 1 ,n 2 ) (25) 



in the continuum limit 



rii 



xi, n 2 



dP(xi,x 2 ,t) 



if 



i) 1 



dt \ dx\ 

Making the change of variable 



dx\ 



dx\dx 2 



P. 



Xi 

1 



x 2 



y 



equation 



V3 
6|) becomes 

dP(x,y,t) 
dt 



(xi 



x 2 ) 



\ dx 2 dy 2 



P. 



(26) 

(27) 
(28) 

(29) 



When either of the domain walls reaches the test site, the 
test site flips. Therefore, if we want P to represent the 
conditional probability that the site has not flipped, we 
must solve equation ( |2q ) with boundary condition P = 
along x\ = and x 2 = 0. This corresponds to solving the 
diffusion equation ( p9| ) with boundary condition P = 
along the lines y = ±a;/v3, for the region x 
x > — 7/v3~ — that is, in a wedge of angle tt/3 with ab- 
sorbing boundaries. 

The diffusion equation is readily solved in a wedge 
of angle ip with absorbing boundaries, and it is known 
that the survival probability for a random walker under 
such conditions decays as t~^/ 2 ^ |14|. After transforming 
equation (p9[) into polar coordinates r and 9, where 



+ y , 

tt/6 



>V2 



tan cj) = y/x, 



(30) 

(31) 
(32) 



we need to solve the diffusion equation for the region 
< 6 < tt/3, with P = on 9 = and 6 = tt/3. The 
appropriate solution, starting from r = r , 9 = 9 is 

3 f°° 

P(r, 9, t) = - V / dX e- xt sin(3m6>) sin(3m6'o) x 

x J 3m (r X 1/2 )J 3m (r\ 1/2 ) (33) 

3 ^ 

= — > sm(3m9) sm(3m9o) x 



where we have used eqn. 6.615 of p5| . To find the 
persistence probability, we need to evaluate Pp{t) — 

r dr JJ 1 ^ 3 d9 P. Performing the integrals (l5| , and tak- 
ing the limit t — ► 00 (where the dominant contribution 
comes from the term m — 1 in the sum) we find 



1 / r 2 X 3/2 



(35) 



On substituting for the initial conditions in terms of 
x = n\ + n 2 , y — (n° — n 2 )/y/3, equation (^5| ) reduces to 
Pp{t) = j^nlnl{nl + nl)t-W. 

In order to calculate the persistence probability, we 
need to average over the possible values of n\ and n 2 . 
The probability in the initial state of a given 'up' spin 
having precisely (rij — 1) consecutive 'up' neighbours to 
its left is (1 — /x)/i™i _1 , so the probability of an initial 
configuration (n^n!]) is (1 — /i) 2 ^™ 1+ " 2 ~ 2 . Summing over 
n\ > 1 and n 2 > 1 gives the final result, 



P P (t) = 



(36) 



One might attempt to perform a similar calculation for 
a site within a majority domain. Here, however, it is im- 
portant to take account not only of the fact that the walls 
of the domain can move, but also that the neighbouring 
domains can vanish before the test spin has flipped. Con- 
tributions from domains that have coalesced will there- 
fore be important, and the necessity to include informa- 
tion about domain-domain correlations makes the calcu- 
lation extremely complicated, if not intractable. 



B. Simulation results 



The persistence probability for both majority and mi- 
nority spins was measured for a range of values of fx. 
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While the measured values of the exponent a for the mi- 
nority domains was found to be close to 1.5 for [i small, 
for fi = 0.5 there appeared to be some evidence that 
the asymptotic behaviour was governed by a different 
exponent. However, the slow onset of the asymptotic 
regime for intermediate values of /i, together with the 
large value of the exponent (> 1, compared with < 1 for 
the Glauber case |13|), led to unavoidably poor statistics 
in the asymptotic regime, and hence it was not possible 
to establish reliable values for the exponent. 
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FIG. 6. Bias plot for the persistence probability of minority 
spins. 

Figure ^ shows a plot of Pp(t)t 3 ^ 2 against time for 
minority spins, for fj, = 0.01. For comparison, the con- 
stant (« 0.2936) predicted by equation (|3g) is also shown. 
There is reasonable agreement of the simulations with 
the independent-domain prediction, suggesting that it is 
a good approximation for small /i, for the time regime 
measured. We interpret the deviation from the constant 
at longer times as more likely to be due to statistics than 
a true systematic effect. Similar data for fi = 0.1 is 
shown, and deviations from the constant predicted by 
(pq) (~ 0.4257) are already quite marked. 
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FIG. 7. Effective exponent for the persistence probability 
of minority spins. 



The effective persistence exponent d(\n Pp) / d(lnt) for 
minority spins is plotted as a function of 1/ In t in Figure 
0, for several values of [i. The purpose of this plot is to in- 
vestigate whether there are any underlying trends in the 
exponent. When a function of the form f(t) oc t a (\nt) b is 
plotted on a graph of this kind, the curve approaches the 
/-axis linearly with gradient b and intercept a. The be- 
haviour for data that converge asymptotically to a power- 
law is for the slope of such a plot to level off to zero in 
the asymptotic regime. For small values of fi, the ex- 
ponent seems to have settled down to a value close to 
— 1.5. However, for /i close to 0.5 there still appears to 
be some trend in the effective exponent, suggesting that 
the asymptotic regime has not been reached. Only for 
\i = 0.5 is there strong evidence that the data converge 
to an exponent different from —1.5. Improved statistics 
and longer times would be needed to give an unequivocal 
conclusion, but we would not be able to do this using 
our current techniques since each set of data required 
approximately 1 week of CPU time. 

The results for the persistence of the majority species 
are even more problematical, because for fi ^ 0.5 too few 
spins were found to have to flipped for the asymptotic 
regime to have been reached. 

The deviations from the free-domain picture must be 
due to the fact that some of the minority domains coa- 
lesce. It is possible to estimate the number of minority 
domains that coalesce from the rate at which majority 
domains vanish. From equation(H), using (pH§), we find 



dN+ 
~dT 



(1 



dl- 



(37) 



z_=o 



where N + = N/2 is the number of minority ('up') do- 
mains per site. Integrating, we find 



dt 



dP- 



dl- 



z_=o 



(i - nY + fx- 



:N+(t). 



(38) 



The term on the left-hand side is the total number of 
majority-domains that vanish between time t and infin- 
ity, which is equal to the number of minority-domain co- 
alescences. Some domains may coalesce more than once, 
so this is an upper bound on the number of minority 
domains that coalesce at least once after time t. The 
fraction of minority domains that coalesce after time t 
is therefore less than /i 2 / [/s 2 + (1 — /i) 2 ] , which vanishes 
like /j, 2 for small [i. 

Measurements from the numerical simulations con- 
firmed the picture that the fraction of domains that co- 
alesce is small for [i small. Nevertheless, it is found that 
the dominant contribution to persistence in the long-time 
limit comes from spins in domains that have coalesced, 
thus explaining the deviations from the free-domain ap- 
proximation. 
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V. SUMMARY 

The low-temperature coarsening dynamics of a driven 
diffusive system - the 1-D Ising model with a driving 
force E satisfying T <C E -C J - has been studied by 
a combination of analytical and numerical techniques. 
Compelling evidence for a mean domain size growing as 
i 1 / 2 , and a domain-size distribution of the form (Q), has 
been presented. These results are exact in the limit where 
one phase occupies a vanishingly small volume fraction 
(i. In the same limit, the persistence exponent for the 
minority phase is 9 = 3/2. The limit of zero volume 



fraction was studied by Lifshitz and Slyozov 16 1 to pre- 
dict the growth exponent (= 1/3) for the case without 
driving force in general dimension, and it is hoped that 
the approach of the present paper might also be usefully 
extended to higher dimensions. 

The random-walk character of the domain dynamics 
suggests i 1 / 2 growth generally, and this is borne out by 
the simulation results (Figures 1 and 2) . The simulations 
also lend strong support (Figures 3 and 4) to the scaling 
distribution (P for general fi. By contrast, the value of 
9 away from the small-/Lt limit is difficult to determine 
numerically, due to slow convergence to the asymptotic 
regime (Figure 7). For [i = 0.5, however, the results seem 
to be inconsistent with 9 = 3/2, suggesting that the ex- 
ponent 9 may depend continuously on /x, as is the case 
for the 1-D Ising model with Glauber Dynamics [ fl2"|| . 

It is tempting to propose an experimental setup 
that might be described by the present one-dimensional 
model. If two immiscible fluids of differing density are 
stirred and placed in a vertical tube, then they will typi- 
cally be able to slide past each other and separate hydro- 
dynamically. However, if the tube is sufficiently narrow 
then it is possible for a state where the denser fluid is 
above the lighter fluid to be metastable; the loss in grav- 
itational potential energy if the interface is tilted slightly 
can be stabilised by the increase in interfacial energy as 
the area of the interface increases. A state consisting of 
alternating quasi-one-dimensional 'domains' can there- 
fore be metastable. There now arises the question of 
adding a noise source to create droplets at the surface. 
Simply shaking the tube will tend to induce tilting of the 
interface, which will lead to hydrodynamic instability. 
The driving vibrations therefore have to be of a wave- 
length much smaller than the width of the tube in order 
to create droplets without exciting the 'sloshing' mode. 
We leave the problems of finding an ultrasonic source of 
high enough intensity, and of preparing the intial condi- 
tion, as challenges for the keen experimenter. 

In a separate paper |l7| we will present results for 
the T = dynamics of a deterministic 1-D scalar field 
model, defined by the modified Cahn-Hilliard equation 



d t( f> = -d 2 x (dl(t> + <f> - (f> 3 ) + E(j)d x cj>. In the small-£ limit, 
the coarsening dynamics of this driven diffusive model 
also exhibit a scaling distribution for domain sizes, and 
a t 1 / 2 growth of the mean domain size. Both models, the 
stochastic Ising model considered here, and the deter- 
ministic model are of great interest in higher dimensions. 
The approach of looking at the limit of small volume 
fraction fx, which proved so successful here, may well be 
fruitful in elucidating the behaviour of these models in 
general dimension D. 
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